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Shegelski, Kavka, and Hynbidai have shown how to calculate the lifetime of a particle 
initially localized in a potential well exactly quantum mechanically by employing a heuristic 
expression with a simple interpretation for the lifetime. Their method allows for the inclusion 
of a bound state, and their results for tunneling through a centrifugal barrier demonstrate 
the major role that bound states play in determining the lifetime. 

However, a drawback to their procedure, especially from the point of view of undergrad- 
uate student accessibility, is that it leads to numerical integrations over time that these 
authors acknowledge are very challenging due to rapid oscillations in the integrands and the 
need to introduce an artificial cutoff for their control. These difficulties are inherent to their 
approach and are not ameliorated by considering simpler barriers such as a square barrier 
or a repulsive 5-function barrier. 

In this note we show that by choosing a slightly different, but equally appropriate, heuris- 
tic expression for the lifetime we can greatly simplify the analysis and arrive at analytical 
results requiring only a single final numerical integration that can be performed quickly and 
without difficulty. We obtain results for the centrifugal barrier which are exact for the case 
where the system has no bound states. For the case where a bound state is present we make 
one approximation which, as our results demonstrate, is relatively mild and preserves very 
well the essential feature of the exact calculation of Shegelski, et al., viz., the influence that 
the presence of a bound state has on the lifetime. 

The expression that we choose to characterize the lifetime is < t 2 > l l 2 as defined by 



and the step function is defined by 9(x) = for x < and 8(x) = 1 for x > 0. 

The wave function at time t is a linear superposition of bound and unbound eigenstates 
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where, as in Ref. (1), 



AP m (t) = [P in {t) - P in (oo)} 6([P in (t) - Pin(oo)]), 



(3) 
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of the Hamiltonian, 



t) = CMx) e- iE »V R + CkM%) e- lEkt '\ (4) 



k 



= V b (x,t) + V u (x,t), (5) 

where, as in Ref. 1, we have assumed the presence of only one bound state. 
Thus, 
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P in {t)= I | m b (x,t)\ 2 dx+ / \^ u (x,t)\ 2 dx + 2Re %{x, t)^ u (x, t) dx. (6) 
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The first term on the right of Eq. ([6]) is simply P in (oo) so that 

pa pa 

Pin{t) - Pm(oo) = / | t) | 2 dx + 2i?e / t)^ u (x, t) dx. (7) 
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The contribution from the second term on the right of Eq. ([7]), the interference term between 

the bound and unbound states,we now argue is relatively unimportant. To see this we need 

only consider the orthogonality property of the bound and unbound energy eigenfunctions, 

L pa pL 

<fil{x)4>k{x) dx = I tpl{x)4>k{x) dx + / (f>l(x)(f)k(x) dx = 0. (8) 
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As illustrated in Fig. 1, the region to the right of x = a is a classically forbidden region 
for the bound state, and its wave function decreases to zero exponentially there. Thus the 
integral over the classically forbidden region (a < x < L) is expected to be quite small, and 
the integral over the allowed region (0 < x < a) must, in turn, be quite small also. On this 
basis we shall neglect the interference term and henceforth represent AP in (t) by 

AP in (t) = [ | V u (x,t)\ 2 , (9) 
Jo 

having noted that in this approximation Pi n (t) — Pj n (oo) > so that the step function in 
Eq. (jSl) is no longer needed. 

We now proceed to the calculation of < t 2 >. From Eq. (10) of Ref. 1 we have^ 

poo 

V u {x,t)= dk<f)(k)sm(qx)e~ i(ka)2(t/t0 \ (10) 
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with 
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and q = q{k) = \Jk 2 — 2Vq (Vq < for a well) as L — > oo and < x < a. The function 
f 2 (k) is defined by Eq. (18) of Ref. 1. For convenience, we record here its actual form for 
a centrifugal barrier with £ = 1: 

f{k) = [(1 + n 2 )a 2 cos 2 a + (1 - k 2 + k 4 ) sin 2 a + a sin(2a)] , (12) 

ft CI 

with k = ka and a = ga. Then, beginning with the x-integration, 

poo poo pa 

AP in (t) = I dk dk'(f){k)(f>{k')e i[{k ' a)2 - (kam/to) / dxsin(gx)sin(g'x), (13) 
Jo Jo 

i /»oo 

dk / dk'®(k, k') e i/(fe ' fe,)t , (14) 



where 



$(MO=0W0(^)xM'), (15) 



,. o' cos(o'a) sin(oa) — 0003(00) sm(q'a) , . 

x(ff, ?o = - — (i6) 



and 



f(k,k')= [(k'a) 2 - (ka) 2 } /t . (17) 

Consider first the numerator in Eq. ([T]). We observe that <&(k,k') is invariant under 
k ^ k' and is an even function of both k and These symmetries allow us to write 
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dtt 2 AP in (t) = \ I dtt 2 AP m (t) (18a) 

II ^ J-OO 
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- dk dk'$(k,k') / dtt 2 e if{k ' k,)t (18b) 
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where the Dirac 5-function was introduced in Eq. (118dj) through its Fourier integral repre 



sentation. Integrating by parts with respect to k and k', and noting that the integrated 
terms vanish,we obtain 

P , 24n , v t n f°° dk f°° „ ,<9 2 ^(fc,fc') rr/J , A r/J , /N1 
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where W(k,k') = $(k,k')/kk'. 
Finally, then 



/*oo rc 

/ dtt 2 AP in (t) = n 
Jo Jo 



dk d 2 ^(k,k') 



(20) 

k'=k 



k dkdk' 

where we have made use of the definition i — 2a 2 in arriving at the last line. A similar 
calculation for the denominator in Eq. ([T]), but without the need to do any integrations by 
parts, results in 

/•CO /-OO 77 

J dtAP in (t)=n J T $(k,k). (21) 

To evaluate < t 2 >, the remaining integration in Eqs. (|20|) and (I2ip must be performed 
numerically, but in each case only a single integration is required. It can be done quickly 
and without difficulty using Mathematica or Maple. 

To compare our results with the exact results of Ref. 1, we recall that for a system that 
decays exponentially, < t 2 > 1 / 2 = ^f2 < t >. Although the behavior of AP in {t) is more 
complicated than a simple exponential, its overall effect on the lifetime should be roughly 
"exponential-like" . Thus, if we define the average lifetime < i > for our approach by 

< t 2 > x l 2 

<i>=^~, (22) 

we expect that < t > will differ very little from < t > the average lifetime as defined in 
Ref. 1. This will enable us to make a direct comparison with the results presented in Fig. 
2 of Ref. 1 by making a similar plot using < t >. For plotting purposes we introduce a 
dimensionless lifetime f =< t > /t and plot f as a function of < e > as defined by Eq. (20) 
of Ref. 1. Our results are shown in Fig. 2 along with the WKB lifetime computed from Eq. 
(28) of Ref. 1. 

A careful comparison of these two plots leads one to conclude that there is no perceptible 
difference between the right-hand half of the dotted curves. Since both calculations are exact 
for the regime in which no bound state is present, we conclude that for comparison purposes 
< t > and < t > may be considered to be the same. With regard to the left-hand half of the 
dotted curves, there is some reduction in the values of the < f > of our results and a sharper 
bend upward as one approaches close to the point where the bound state disappears. It is 
in this region that one expects the interference term, which we have neglected, to exert its 
greatest influence. But qualitatively, our results for the bound state regime are remarkably 
similar to the exact results of Ref. 1. Thus, they too represent a significant improvement 
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on the WKB approximation while, at the same time, maintaining a much greater simplicity 
of computation. 

Finally, for those who wish to avoid the use of the Dirac ^-function, we offer an alternate 
method for calculating the numerator and denominator in Eq. ([T]). It consists of inserting 
a convergence factor e~ at into the original expressions for these quantities so that the t- 
integration can be carried out immediately (along with the x-integration). The k '-integration 
is then performed using the method of residues after which one can safely set a —>■ 0. There 
then remains only the k- integration that, once again, must be done numerically. 

Problem 1: Apply the convergence factor /residue method to the calculation of the de- 
nominator in Eq. ([T]) and show that the result obtained is identical to that given by Eq. 

Problem 2: Apply the convergence factor/residue method to the calculation of the numer- 
ator in Eq. ([1]) and show that the result obtained is equivalent to Eq. (|20|) by demonstrating 
that it yields exactly the same numerical result as Eq. ( l20l) . 



* Electronic address: intemann@temple.edu 



1 Mark R. A. Shegelski, Jeremy Kavka, and Jeff Hnybida, "Quantum mechanical versus semiclas- 
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2 We have set m = fi = 1 at this point for convenience. 
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FIG. 1: The potential for the centrifugal barrier problem. The dashed curve illustrates a typical 
bound state wave function for the case where a bound state is present. 
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FIG. 2: The tunneling time as a function of the dimensionless energy < e >=< E > /V&, where 
< E > is the expectation value of the particle's energy and Vf, is the height of the centrifugal 
barrier. The solid circles are the results of the present work and the continuous curve represents 
the WKB approximation. 
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